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Abstract. We assume that the Universe has a non trivial topology whose compact spatial sections have a volume 
significantly smaller than the horizon volume. By a topological lens effect, such a "folded" space configuration generates 
| multiple images of cosmic sources, e.g. clusters of galaxies. We present a simple and powerful method to unveil non 

■ ambiguous observational effects, independently of the sign of the curvature and of the topological type. By analogy 
with techniques used in crystallography, we look for spikes in the pair separation histogram between cosmic objects 

■ in 3-D space. The spikes due to multi-connectedness should stand out dramatically. Moreover, their positions and 
I | their relative amplitudes would be definite signatures of the topological type and of the underlying geometry. Such a 

^» ■ statistical method would thus reveal the shape of space. As illustrative examples, we perform numerical simulations in 
| fl = 1 Friedmann universes with the six possible closed orientable topologies, which prove the efficiency of our method. 
*^ Presently available 3D catalogs of galaxy clusters are not deep enough to test our method at sizes greater than lower 
| limits w 600 hr 1 Mpc previously obtained by other methods. With extensive redshift surveys currently in progress the 

■ situation may change in the next decade. 

so : 

0^ ■ Key words: cosmology: large scale structure of Universe, topology 

"o 

CF; 

• • ■ 1. Introduction 

> ■ 

J^j ' It is presently believed that our Universe is correctly described by one of the spatially homogeneous and isotropic 
Friedmann-Lemaitre models. Their spatial sections are of the elliptical, Euclidean or hyperbolic type according to the 
C$ ' sign of their constant spatial curvature. Most studies in the field assume that the topology of space is simply-connected, 
but the possibility of a non trivial topology becomes an increasingly popular topic in theoretical and observational 
cosmology (for a review, see Lachieze-Rey and Luminet 1995, hereafter LaLu95 and references therein; also Fagundcs 
1995). 

For the simply-connected models, the finite or non finite character of space is linked to the sign of spatial curvature: 
elliptical models have finite volumes whereas Euclidean and hyperbolic models have infinite volumes. In that case the 
question of the extension of space (one of the oldest cosmological problems, going back to almost twenty five centuries) 
is reduced to the estimation of the average mass-energy density of the Universe (including a possible cosmological 
constant's contribution). The situation is different for multi-connected models (Ellis 1971), since the three geometries 
of constant curvature admit compact space forms (although the elliptical family do not admit non compact ones). 

Present observations do not clearly indicate the value of the spatial curvature of the Universe, neither whether 
space is simply or multi-connected. However observational effects are expected if the physical space is compact in 
at least one direction and if its corresponding size is smaller than the horizon distance. We call such a space "small 
folded" . Two of us (LaLu95) published a critical review of the various methods already proposed to check the cosmic 
topology; most of them present at least one of the following drawbacks: (i) they apply to one particular type of multi- 
connectedness only, for instance the torus like topology; (ii) they rest on strong assumptions about the cosmological 
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model, for instance an Einstein-de Sitter universe, or on the properties of a peculiar population of cosmic objects, for 
instance quasars; (iii) they do not provide an unambiguous signature of multiconnectedness. 

In this article we propose to test if space is a small folded one by studying the histogram of pair separations 
between cosmic sources. Our method is an improved version of the construction of the correlation function of clusters 
already pioneered by Schvartsman and his group (references in LeLa95). It is free from the preceding drawbacks and 
is independent of preliminary assumptions about the topological type. 

We recall general properties of multi connected universe models in section 2. We present the pair separation 
histogram method in section 3. We check its validity on simulated universe models in section 4. We perform the test 
on a 3D catalog of galaxy clusters and draw conclusions in section 5. 

2. Small folded universes 

2.1. General properties 

Mathematically, a space is simply-connected if every loop can be continuously shrunk to a point. If not, the space is 
multi-connected. A multi-connected space is conveniently described by its fundamental polyhedron F and its holonomy 
group r. F is convex, with a finite number of faces which are identified by pairs. T is generated by transformations 
which carry a face to its homologous one. The latter are isometries without fixed point. The fundamental polyhedron 
is transformed into its images 7F by the holonomics 7 G T. The reunion of the 7F's form a regular tiling of the 
so-called universal covering space. A simply-connected space is identical to its universal covering space, whereas 
a multi-connected space is the quotient of the universal covering space by the holonomy group of its fundamental 
polyhedron. 

Applied to the spatial sections of cosmological models, the fundamental polyhedron indeed identifies with the 
physical space, in which objects like galaxies or clusters are located, whereas the universal covering space identifies 
to S 3 , IR 3 or H 3 according to the sign of the curvature (we assume, as usual in cosmology, homogeneity and local 
isotropy of space). 

A given cosmic object lies at a given position in physical space, i.e. in F. Its images under T lie in the universal 
covering space. The latter is thus the "observer's world", which may drastically differ from the real world. As an 
example, the observer may see several images of the same cosmic object if F is smaller than the particle horizon 
distance in at least one space direction. This gives rise to the appellation of "small folded" universe. By convention 
we call the nearest source the "original" and the other images "ghosts" . 

Small folded universe models are easier to understand with the following remarks: (i) although it is not the 
physical space, the universal covering space in a given multi-connected Friedmann-Lemaitre model has exactly the 
same properties than the physical space of the corresponding simply-connected model; (ii) a ghost object in the 
universal covering space of a multi-connected model has the same properties (observed distance, rcdshift, age) as the 
object located at the same position in the real space of the corresponding simply-connected model. 

A multi-connected universe model is characterized by some spatial scales associated to the fundamental polyhedron. 
Let us call a its smallest length. In a non flat space, characterized by its present curvature radius i?o, the ratio a/Ro is 
geometrically constrained to specified values (LaLu95). For instance it has a maximum (resp. minimum) value if k > 
(resp. k < 0), whereas it remains arbitrary in flat space. The fundamental polyhedron also involves another scale (3, 
the maximum length inscriptiblc in it (this is for instance the diagonal for a parallepipedic fundamental polyhedron, 
which characterizes the hypertorus) . This is also the maximum distance between 2 images of the same object belonging 
to adjacent cells. 

Directly observable effects are expected if a or (3 are smaller than the horizon size. By construction, only original 
images (no ghosts) are present up to a distance a/2. Between a/2 and (3/2 there is a mixture of original and ghost 
images. Beyond (3/2 there are only ghosts. Of course, the smaller the basic cell, the easier to detect topological effects. 
The number of potential ghosts depends on the cosmic parameters and on the topology. For a given object, there are 
as much potential ghosts as cells in the universal covering space; their number is thus finite in the case of positive 
constant spatial curvature, infinite otherwise. An obvious limitation is that only ghosts nearer than the particle horizon 
may be seen, so that their number is at most equal to the number of cells within the particle's horizon. In addition, 
the number may be drastically reduced by the fact that in practice, we can only observe below some magnitude or 
rcdshift cut-off, depending on the type of objects (galaxies, clusters, etc.), on the instrumentation and various other 
limitations. 
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2.2. Observational constraints 

No direct observation presently indicates that our universe is a small folded one. By direct observation we mean 
recognition of ghost images of known cosmic objects or configurations. This places lower limits on the size of the 
physical space according to the type of objects. Some tests involve the diffuse cosmic microwave background, other 
the discrete sources such as quasars or clusters. From the COBE/DMR results, Stevens et al. (1993), Starobinsky 
(1993) and de Oliveira-Costa & Smoot (1995) have claimed to rule out multi-connectedness on sub-horizon scales. 
However their results rely onto various disputable assumptions, for instance the absence of strong reionization after 
recombination (see section 12.3 in LaLu95 for an extended discussion). Fagundes (1995) tried to fit the cold and hot 
spots in COBE/DMR maps of the cosmic microwave background with the predictions of some compact hyperbolic 
models with £1 < 1, whereas Jing and Fang (1994) fitted the two-point angular correlation function of the cosmic 
microwave background. At present time, limits derived from the cosmic microwave background observations are not 
model independent and cannot be considered as definitive yet. Thus it remains valuable to check the cosmic topology 
with a different method. 

Tests which involve discrete sources must consider populations of objects extending deep enough in space to check 
large dimensions. Quasars seem potentially interesting, but their estimated lifetime is short compared to the expected 
time necessary for a light ray to turn around a small folded universe: a typical lifetime less than 10 8 yrs would only 
allow to investigate scales smaller than 200 Mpc. 

The observations of clusters and supcrclusters of galaxies seem more adapted. They have led to the limits a > 
60 h- 1 Mpc (Gott 1980) and (3 > 600 h' 1 Mpc (Sokoloff and Schvartsman 1974, Fang & Liu 1988), where h is the 
Hubble constant in units of 100 km sec -1 Mpc -1 . To our knowledge, they are the best observational constraints 
derived from discrete sources. They leave room for many observable effects in a folded universe. To get a qualitative 
idea, let us assume an Einstein-de Sitter model (zero curvature) with a cubic hypertorus of length L as fundamental 
cell. The redshift-distance relation is d(z) = du (l — (1 + z) -1 / 2 ), where d(z) is the distance in the universal covering 
space, and dn = 2c/ Hq 6000 h^ 1 Mpc is the horizon distance. Table 1 gives the multiplication factor, namely the 
number of ghost images of a given original object, for various cell sizes compatible with present observations. 



Table 1. Multiplication of images in a cubic hypertorus universe with fundamental lengths L, in units of h' 1 Mpc. z a is the 
redshift until which there are only originals, z@ the redshift beyond which there are only ghosts, M(dfr) the multiplication factor 
until the horizon distance du, M(z) the multiplication factors until various redshift cut-offs. Magnitude limitations, absorption, 
luminosity evolution and other observational biases are neglected. 
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M(l) 


M(0.5) 


500 


0.09 


0.16 
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1200 


180 


45 
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0.19 


0.37 
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150 


23 
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0.31 


0.63 


279 


45 
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1.5 


2000 


0.44 


0.98 


110 


20 


3 




2500 


0.60 


1.45 


60 


10 


1.5 





The limits mentionned above are mainly derived from the absence of ghosts of peculiar objects (the Milky Way, 
some known clusters or superclusters). However some ghosts would be unobservable for various reasons not linked to 
the geometry (absorption by dust, difficulty to recognize the object, etc.). This motivates a search for statistical tests 
which would be less dependent on peculiar objects and would exploit the largest amount of the information contained 
in 3-D catalogs. 

Different methods have been proposed, such as the search for a periodic or quantized distribution of objects in 
redshift or in distance (Fang 1990). For instance, it may be tempting to invoke multi-connectedness for explaining the 
apparent periodicity observed by Broadhurst et al. (1990) in a pencil-beam catalog of galaxies. But Park and Gott 
(1991), among others, were able to reproduce comparable results in ordinary simply-connected universe models. As 
emphasized in LaLu95, such tests cannot be decisive, and our numerical simulations below confirm this point of view. 

3. A pair separations histogram for testing cosmic topology 

Since the above methods do not provide a clear signature of multi connectedness, we propose a more general one, 
independent on the peculiar topology or on the peculiar cosmic objects. In this purpose, we go back to the basic 
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property of folded universes: each image of a given object is linked to each other one by the holonomies of space. 
Without knowing the geometry or the topology we do not know what are these holomies, but we know that they are 
isometries. Therefore, any pair of images comprising an original and one of its ghosts, or two ghosts of the same object, 
reflects an isometry of space. We call such a pair a gg-pair. A signature of multi-connectedness would be that, among 
all pairs of images, a significant proportion are gg-pairs. So the question is: how to extract gg-pairs from ordinary 
pairs ? The answer is provided by the histogram of space separations between all pairs of images in a 3-D catalog, as 
we see now in more detail. 

A topology is characterized by its holonomies, which are combinations of the generators {"fk} of the holonomy 
group. Each generator -fk is itself an holonomy, to which is associated an identification length A& related to the size 
of the fundamental polyhedron. For instance, in toroidal models the A& are equal to the edges L e , L a and L u of 
the parallepipedic fundamental polyhedron. More generally, the A| are related to L\, L\ and L\ by linear relations 
involving integer coefficients. 

All the gg-pairs corresponding to the same generating holonomy -fk are characterized by the same separation A^ 
(in proper comoving distance units). Thus, if we draw the histogram of the pair separations, or rather the squared 
separations, these gg-pairs will emerge from ordinary pairs as a spike located at the position A|. In addition to the 
generating holonomies, other holonomies appear as compositions of them, also giving characteristic spikes. They also 
generate gg-pairs with squared separations 



where the Nk are integers. 

The principle of our test is to recognize the presence of spikes associated to these values, in a catalog of observed 
cosmic sources. Let us assume that our catalog of N objects (thus containing N(N — l)/2 pairs) has a characteristic 
volume about F times that of the fundamental polyhedron. The latter thus contains about N/F (original) objects, the 
other being ghosts. Each ghost gi (excepted near the edges) is transformed by the holonomy generator -fk to give an 
other ghost 32, their separation being A&. Thus, we expect about N (in fact less because of edge effects) corresponding 
gg-pairs with separation Afc, which will emerge as a spike above the background contribution of ordinary pairs. If 2 
or 3 identification lengths are equal, this number is to be multiplied by 2 or 3 and the corresponding spike will be 
enhanced. 

In a small folded universe, the spikes would be observable and provide a signature of multi-connectedness, like in 
a crystallographic lattice. The integer values Nk would be measurable and (together with the relative heights of the 
spikes) indicate exactly the type of topology involved. This is illustrated in the next section. Note that the practical 
calculus of separations depending on the curvature of space, the analysis has in principle to be performed in the three 
different cases k = — 1, 0, 1. However, the Euclidean approximation applies well when the catalog does not extend too 
deeply. 

4. Simulated folded universes 

Locally Euclidean models 

In order to check the validity of our test, we have generated simulated catalogs of a distribution of objects (typically, 
galaxy clusters or quasars) in a multi-connected universe. For simplicity, we chose a flat Friedmann-Lemaitre (Einstein- 
de Sitter) universe, whose universal covering space is 1R 3 . We investigated the six compact oriented topologies, denoted 
71 to 76 (see Wolf 1984, and Fig. 17 of LaLu95 for a pictorial representation). Each type is specified by a fundamental 
polyhedron with three edge lengths L e , L a , L u , and the structure of its holonomy group. The latter may be characterized 
by the formulae identifying the coordinates r = (x, y, z) and r' = (x' , y', z') of a gg-pair : 
for types Ty, T 2 , T 3 




(1) 




(2) 



for types T4 




(3) 
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for types T 5 , T e 

(I -1/2 0' 

r' = RoMr) r + V3/2 | ( n a L a ) (4) 
\0 1 

where Rox(9), Roy(&) and Roz{9) are the rotation matrices of angle # about the three coordinates axis, n e , n a and 
n u are integers denoting the position of a cell in the universal covering space, L e , L a , L u and on are given in Table 2. 




Table 2. The six locally Euclidean, closed, oriented 3-spaces 



Type % 


Fundamental polyhedron 


OLi 


Tx 


parallelipiped, L e , L a , L u 





T 2 


parallelipiped, L e , L a , L u 


Tvn u 


% 


parallelipiped, L e = L a , L u 


tt/2 n u 


% 


parallelipiped, L e , L a , L u 




% 


hexagonal prism, L e = L a , L u 


27r/3 n„[mod 3] 


% 


hexagonal prism, L e = L a , L u 


27r/6n„[mod6] 



Type Ti corresponds to the well known hypertorus. Types T3, T§ and Tq constrain two lengths to be equal. The 
volume of a cell (also of the physical space) is V — L e L a L u for types T1-T4, and V = 3V3/2 L 2 e L u for types T§ and 
%. 

4-2. Analyses of simulated catalogs 

We randomly distribute 50 cosmic objects within the fundamental cell and we calculate the ghost images in the 
universal covering space, to simulate the appearance of the sky up to a redshift cut-off z = 4. For galaxy clusters, we 
neglect the peculiar motions which could slightly deviate the positions of ghosts from a rigorous distribution. 

In a first series of simulations we have chosen equal edges L e = L a = L u , adjusted so that the fundamental cell 
has a volume (1500 Yr 1 Mpc) 3 . According to Table 1, the number of ghosts is about 45 times the number of originals. 
Ghosts appear at z > 0.31 and originals disappear at z > 0.63. Figure 1 shows the 2-dimensional appearance of the 
sky. Fake large-scale structures are generated but, as expected, no information about multi-connectedness appears. 

We have also plotted in Fig. 2 the histogram of redshifts, calculated for the hypertorus model. No periodic pattern 
nor any sign of multi-connectedness appears. This confirms former simulations by Ellis and Schreiber (1986) and the 
statements by LaLu95 according to which such histograms cannot provide a valuable test for cosmic topology. 

Next, for each simulated catalog of objects we have drawn the pair separation histogram. They are displayed in 
Fig. 3, from which the information about multi-connectedness and the topological type spectacularly springs out. 

Spikes, which gather all gg-pairs linked by holonomies, appear at squared separations Aj obeying to Eq. (jl|). A 
straightforward calculation of the separations between (x',y',z') and (x,y,z) according to Eqs. (2-4) and Table 2 
shows that Eq. (1) may be rewritten under the very simple form 

A 2 

n 2 e +n 2 a + nl = ^73 (5) 

for the triplets on integers (n e , n a ,n u ) compatible with the topology. The values of n span a range limited by the finite 
size of the catalog (0 to 5 in our calculation). For a given topology the positions of spikes are displayed in Table 3. 

For a better understanding of Table 3, let us comment in more details the spectrum of spikes obtained for types 
T\ (cubic hypertorus) and T4. For 71, the holonomies are generated by the three translations along the edges of 
the parallepipedic cell, with magnitude L — V 1 I Z . In such a simple lattice structure, all the triplets of integers 
(n e ,n a , n u ) e {0, . . . , 5} realize an holonomy. Spikes thus appear at all locations Aj such that A^/V 2 / 3 is the sum of 
three squared integers. Most integers except 7, 15, 16. . . are solutions. For T4, the holonomies are generated by two 
translations and a glide reflection. For such a lattice structure, only a few triplets of integers (n e ,n a ,n u ) realize an 
holonomy, so that Eq. (||) has much less solutions. 

Let us also examine the amplitude of the spikes. For type Ti for instance, the spike at position Af/V 2 ^ 3 = 1 is 
generated by the contributions of the triplets (1, 0, 0), (0, 1, 0) and (0, 0, 1), whereas the spike at position A 2 /V 2 ^ 3 = 5 
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Fig. 1. Appearance of the sky (equal area projection) in an Einstein-de Sitter universe with topology T\. The fundamental 
polyhedron is a cubic hypertorus whose size is 1500 h^ 1 Mpc. 

is generated by the triplets (2,1,0), (2,0,1), (1,2,0), (1,0,2), (0,1,2) and (0,2,1). Thus, spike 5 is more intense than 
spike 1. All the characteristics of Table 3 can be described in a similar way. We conclude that the very existence 
of spikes reveals the multi-connected nature of topology, whereas their positions and amplitudes discriminate the 
topological types. 

We also performed the simulation in a toroidal model like 7i , but with unequal identification lengths. Figure 4 shows 
the result for L a , L e = V^L a , and L u — \/3L a , so that V — (1500 h^ 1 Mpc) 3 . The lengths being not commensurable, 
the spikes are more numerous, corresponding to all positions fundamental lengths and there multiples, but less intense 
since no gg-pair corresponding to different A& accumulate in a same spike. 

Our simulated catalogs are however highly idealized so that the spikes due to gg-pairs stand out dramatically above 
the background distribution of other pairs. In a more realistic distribution, it would be more difficult to recognize the 
spikes, since various effects would contribute to spoil the sharpness of the spikes and decrease the signal to noise ratio. 
For instance, the real catalogs do not include the region masked by the galactic plane, and cover only a small area of 
the sky. In order to examine the consequences of this, we also created a catalog, in a simulated small folded universe, 
which is of limited solid angle. We observe that the signal fades out when the aperture angle of the catalog goes down 
to about 20° (see Fig. 5). 

Also, clustering in a simply-connected universe can generate noise spikes that mimic spikes due to multi- 
connectedness. For instance, the spikes in the pair separation histograms shown in the Broadhurst et al. sample 
(1990) and in the Parks and Gott (1991) simulations are due to noise. More precisely they are caused by the large 
number of pairs relating the many galaxies in a first cluster to those in a second cluster, all at the cluster separation. 
Generally, such spikes currently also occur in N-body simulations that show clustering. 

5. Results and discussion 

5.1. Analysis of a cluster catalog 

We have applied our test to a 3-dimensional catalog of galaxy clusters compiled by Bury (private communication). 
This catalog includes all Abell and ACO clusters with published redshifts. It contains 901 clusters, up to a maximal 
rcdshift z max « 0.35, corresponding to 840 /i -1 Mpc in an Einstcin-de Sitter universe (although only 12 objects have 
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Fig. 2. The redshift histogram generated by 50 randomly distributed objects in the fundamental cell of a 71 universe. 



Table 3. The spike's spectrum of small folded Eucidean universes. The rank of the three stronger spikes are denoted by the 
exponant. 
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Fig. 3. Histogram of squared separation distances between all pairs of images for the 6 topological types, with equal lengths. 
The spikes reveal repetition scales related to the size of the fundamental polyhedron. 
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Fig. 4. Histogram of squared separation distances for a Ti universe with unequal lengths. 



Fig. 5. Histogram of squared separation distances inside a conical catalog of aperture angle 20°, for a Ti universe with equal 
lengths. 



z > 0.26). Due to obscuration by the galactic plane, the shape is a double cone of aperture ~ 120°, sufficient to allow 
the desired effects to be observable. Figure 6 shows its projection onto the sky. 

Figure 7 presents the histogram of separations where two suspicious spikes are detected corresponding to separation 
of 270 ft. -1 Mpc and 382 hr x Mpc. The ratio of these quantities is very near \/2. If those spikes are really due to 
topology, this means that we are seeing multiple images of some given objects. In order to check this hypothesis, we 
have first simulated the Bury catalog: we distributed ramdomly 30 objects in a fundamental cell with L e = L a = L u = 
270 hr 1 Mpc, put a redshift cut-off at z = 0.26, and limited our set of ghosts by a double cone of aperture 120°. The 
number of real objects is chosen to obtain 901 ghosts in the final catalog. Then, we plot the histogram of separation 
on Fig. 8. 

The comparison with Figs. 7 leaves no room to doubt. Multi-connectedness at a scale w 270 /i -1 Mpc would give 
a much stronger signal than observed. To check further, we also selected all pairs contributing to the observed spikes. 
Were they due to multi-connectedness, then the associated separation vectors would clump in preferred directions, 
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Fig. 6. Two dimensional equal area projection onto the sky of the Bury cluster catalog. 



Fig. 7. Histogram of squared separation distances for the Bury cluster catalog. 



corresponding to the principal directions. Clearly, this is not the case as indicated by Fig. 9. Thus the observed spikes 
are only due to noise and no effect of multi-connectedness appears. Given the depth of the catalog, we conclude 
to the limit a > 650 h~ x Mpc, comparable to those already obtained from large scale objects distribution. However 
applications of our test will be decisive when deeper 3D catalogs are available. Current observational programs devoted 
to extended redshift surveys will offer this possibility within the next decades. 

5.2. Conclusion 

We have presented a new method to check and characterize a possible multi connectedness of space on sub-horizon 
scales. Our test is free of specific assumptions about the geometry and the topological type. It is based onto the 
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simple fact that a displacement carrying a ghost image to another of the same object is an isometry in the universal 
covering space. Thus, in a multi-connected space, we expect that the histogram of all pair separations between objects 
of a large 3-dimcnsional catalog exhibits strong spikes, corresponding to combinations of the identification scales. 
Numerical simulations performed on Euclidean "small folded" universe have shown that the relative positions and 
amplitudes of the spikes characterize in an unique way the topological type: different topologies generate different 
recognizable spikes structures. Applying this test to present 3-dimensional catalog of galaxy clusters does not provide 
new limits on the lower size of the universe. 

Taking into account various effects such as luminosity evolution, peculiar velocities and clustering of real objects 
within the fundamental cell would not change the value of our test. In the future, it will be interesting to apply the 
test to more extended catalogs, and to collections of objects of different types. Note also that the test can be made 
still more efficient by considering pairs of specific configurations rather than isolated objects, like strings of images 
with a given shape. 

In the case of a significative distorsion, the pair separation histogram would also provide a powerful and purely 
geometrical method to determine the sign of the curvature of space. This results from the fact that the crystallographic 
structures of the holonomy groups differ for the cases k = 1,0,-1. Thus, beside its own specific interest, multi- 
connectedness could bring an answer to a long-standing problem of observational cosmology. 

Acknowledgements. We thank the anonymous referee for useful comments on the first manuscript and Bury for providing us his 
clusters catalog. 
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